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SUMMARY 


j NASTRAN has been used to solve two types of electromagnetic field problems. 
[The diffusion equation and the boundary conditions valid for problems of these 
kinds together with a replacing potential energy function have been given. The 
I extent to which an analogy with finite element displacement and temperature 
approaches holds is indicated. The outputting of complex quantities is made 
' possible after adjustment of standard rigid format 1 input data blocks to 
, module SDR2. The applications made involve the study of the proximity effect 
in a system of three parallel conductors and the analysis of the magnetic 
' field in the vicinity of the points of contact in circuit breakers. 
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INTRODUCTION 


Obviously the knowledge of electromagnetic field distributions is of prime 
interest in Electric Power Engineering. 

I Transient magnetic fields present in conducting materials cause the introduction 
of eddy currents in such materials. Energy losses leading to temperature in- 
creases result from these currents. The most familiar application of magnetic 
fields is in transformers where they determine the modification of current and 
voltage levels. They have also proven to determine the design of advanced types 
of power generators such as fusion reactors and MHD generators. 

These and many other reasons explain why the analysis of electromagnetic 
, problems is as important as it is. After finite difference techniques were used 
for these problems originally we now observe more and more interest in finite 
element methods. Among the applications of finite element techniques the works 
of Chari and Silvester (Ref.|l|), Chari (Ref.|2|), Silvester and Rafinejad 
(Ref.|3|) and Donea et al (Ref.|4|) have to be mentioned. 


I 

I 

I 

[ 

1 


557 


Rather than writing a stand-alone finite element code for the exclusive 
study of electromagnetic field problems it was felt worthwhile to investigate 
the applicability of the NASTRAN system. ! 

When studying the equations and boundary conditions describing the problems 
involved it shows out there are many similarities with the ones NASTRAN is de- 
signed to deal with. The equation of balance valid here is identical to the j 
equation of heat balance in solids except for the unknown quantity which is a 
vectorial instead of a scalar. When considering harmonic oscillations the tran- i 
sient term in the equation of balance may be transformed to a stationary one and 
a diffusion equation is obtained. In case the unknown field was a scalar quan- 
tity the solution to the magnetic field problem would be similar to the one | 
valid for one-group neutron diffusion. 

A NASTRAN solution to the neutron problem is reported in Ref . | 5 | . i 

Electromagnetic field problems discussed in the present work mean a class of 
other nonstructural applications of NASTRAN and the approach applied here in - 

fact means an extension of the research discussed in this reference. ! 


FORMULATION OF THE ELECTROMAGNETIC FIELD EQUATIONS 


Electromagnetic problems are described by Maxwell’s equations. 

Without going into details these relations will be restated here for complete- 
ness. 

Placing an electrically conducting body in a harmonically oscillating magnetic 
field B will cause eddy currents to be introduced in it. The relation between 
B and the electric field E introduced is given by Faraday’s law of induction: 


V X E 


at 


( 1 ) 


Current density J is obtained from Ohm’s law: 

J = aE (2) 

where a denotes the conductivity. 

These currents generate a magnetic field which interferes with the original 
one. When restricting ourselves to slowly changing systems (f < 100 Hz) we may 
use the relation: 

V X B = y J (3) 

where y stands for the permeability of the conductor. _ 

We introduce the concepts of the magnetic vector potential A which is related 
to B by: 

B = V X A W 

and the one of the electric potential ij; defined by: 
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(5) 



= - Vi|; 


E indicates the electric field present when magnetic influences are omitted, 
el 

by definition 


j V.A = 0 

Moreover 


( 6 ) 


” 3t * 


where represents the current caused by the electric potential. 
Equations (1) ... (5) can be combined to obtain the diffusion equation 


(7) 


v2a = ya - yJ^ 


( 8 ) 


jsince we are dealing with harmonic systems it is favourable to replace A by 
"Ae^^^^ and7g by 7 with j standing for operator /^. This will make 

;Eq. (8) read: 


' V^A - jtoyaA + yJ^ = 0 (9) 

I 

^ which describes a stationary rather than a transient state. 

^ However we have to pay a penalty since we are now dealing with complex quan~ 

. . dk . 

The boundary conditions valid in most cases are reflective (-^ = 0) or consist 

of prescribed values of X. For certain types of problems the location of the 
boundary may not be clearly established and is in fact situated at infinite 
distance from the area of interest. Moreover we have to meet the requirement 
I that eddy currents should not give any contribution to the total current. This 
means the eddy current density integrated over the volume V vanishes, i.e. 

’ Ja dV = 0 (10) 

V 

I 

* The solution of Eq. (9) is obtained by solving the equivalent variational 
I problem, i.e. by imposing the requirement of the relevant energy func- 
tional P to be stationary. 

[ 

I 


P = ///<i A^)dV + ///yj A dV = Stationary 0 0 

When contrasting' the first term of this functional with the thermal potential 
function it shows out to the neglect of a constant factor this term contains 
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the thermal potential function. Moreover when omitting a constant factor the ■ 
capacitive energy function is implied in it. 

This indicates electromagnetic field problems may be computed using finite 
element heat conduction computer programs at least modified in such a sense th^ 
the heat capacity matrix is added to the heat conduction matrix. Alternatively 
finite element structural mechanics codes with provisions to account for aniso- 
tropic material properties may be used. Specification of particular material i 
properties and addition of the mass (or damping) matrix to the stiffness matrix 
is required in enabling the application of this alternative approach. 


THREE PHASE BUS BAR SYSTEM 


As a first example on how NASTRAN can be applied to deal with these prob- 
lems we will discuss the interaction of electromagnetic fields in a system of j 
three parallel conductors carrying alternating currents with mutual phase lags 

T . 27T , 

equal to -y. 

For this type of problem the following assumptions are valid: ! 

i 

1) The magnetic vector potential and the source current densities in the bars 
have only components in the longitudinal direction of the system. They are | 
invariant with this direction and vary sinusoidally with time, so the prob- 
lem is essentially two-dimensional. 

2) The fields are assumed quasi-stationary , so that displacement currents may : 
be neglected. 

3) Hysteresis, magnetic saturation and temperature effects of resistivity are 
negligible. 

Because of the first of these assumptions A and J may be considered to be 
scalars. At infinite distance around the system this quantity vanishes and in 
order to apply this boundary condition properly it seems desirable to make use 
of something like the apparently contradictory concept of infinite finite el- 
ements recently introduced by P. Bettes, University of Wales. The lack of this 
type of elements can be overcome by restricting ourselves to a usual finite mesl 
(Fig.|l|) and diminishing of the results with a certain quantity, resulting froi 
Eq. (10). _ 

The phase lagged volumetric quantity yJ^ present in Eq. (9) is composed of thre< 
sets; one for each of the bars. Since the unknown is a scalar quantity and sinc^ 
in the displacement approach volumetric loading (GRAY bulk data cards) can nei- 
ther be computed to obtain complex quantities nor be varied over subcases it is 
appropriate to use the NASTRAN heat conduction capability where we can use QVOL 
cards. Complex values of volumetric heat generation rate cannot be specified. 
This deficiency can be overcome by stepwise generation of the electromagnetic 
load in three different subcases and subsequent phase lagged superposition of 
the load vectors, making use of complex multiplication factors specified in the 
parameter section of DMAP module ADD. 

The load vectors have to be computed separately instead of being appended. 

This can be accomplished by specification of different single point constraints 
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I 

&r each of the individual subcases. Therefore an SPOINT has been introduced 
|ith three differing prescribed values. 

li order to determine the solution at the right level the subsidiary condition 
Eq. 10) has to be applied. This means all elements of the solution vector {u } 
ave to be subtracted by its mean value: ® 


u 


m 




02 ) 


here S stands for the cross-sectional surface of the conductors. This multipli- 
ation factor is computed from: 


~ {u }'^|m |{I } 

pS g ' gg' g 

|ith {Ig} denoting a vector of order g with all elements equal to !• 


(13) 


!ince SPR's 458 and 483 have not yet been corrected in Level 15.5 of NASTRAN 
^odule ADD is unable to add complex input matrices. Consequently, at present the 
subtraction mentioned has to be performed either by hand or in a stand-alone 
togram. For the same reason the computation of the current density can not be 
lone in NASTRAN. 

ipecial care has to be taken that complex quantities are output from module 
DR2. This requires reformatting of some of its input data blocks. In order to 
Lccomplish this a dummy table of eigenvectors CLAMA has been specified on DTI 
>ulk data cards. A dummy EIGC card was used in order to have IFP create a 
)YNAMICS data block* Moreover modules CASE and DPD were introduced to provide 
CASEXX, EQDYN and SILD required to make SDR2 print complex eigenvectors. For the 
purpose of this application this means the printing of complex nodal vector po- 
:ential values. The rigid format Alter resulting from these modifications is 
;iven in Appendix I. 


The range of parameters used for the computation of A reads as follows: 


specific conductivity 
)ermeability 

ilectrical current density 

^ngular frequency 

bis of total current through each 
[ conductor 

::ross-sectional surface of one 

conductor 


a = 3. 7 10^ 0hm“^ m ^ 
y = 47 Tio”^ Henry m“^ 
Jg = 8.71 1 A m"^ 

0 ) = 314.16 rad. sec” 

J = 1500 A 

S = 1722io"^ 


rhe distribution of the current density J derived from A demonstrates the large 
impact on the uniform distribution due to the presence of the magnetic field, 
rhis phenomenon is called the proximity effect. 

?or cross-sectional surface S and specific electric conductivity a the heat pro- 
iuction per unit length is computed from 
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(U) 


Q = // i 17^1 ds 
s 

Obviously, for any fixed ^et current through the system, the minimum value of Q 
(=0^^^) is obtained when J is uniformly distributed. 

The economy of the conducting system is determined by the ratio of Q over 
which is called the resistance ratio. Q cannot be computed in the DMAP 
sequence since SPR's 458 and 483 prohibit the derivation of J (Eq. 7). A stand- 
alone calculation of Q is simplified when making use of j 

ADJ2 = {Ug}^lMgg|{Ug}. 05) ^ 


referred to in the ALTER packet. 

To the neglect of the heat production Q* can be computed from 



(16) 


with {u 


*jT 


standing for the transposed complex conjugate of 


Q* = /{Re(ADJ2)}2 + {Im(ADJ2)}2 


(17) 


For the present configuration the resistance ratio was found to be 2.4. 


MAGNETIC FIELDS IN CIRCUIT BREAKERS 


One of the aspects to be taken into account in the design of circuit 
breakers is the phenomenon of electrical discharge. 

The occurrence of this phenomenon depends on the amplitude and phase angle of 
the magnetic field present in the neighbourhood of the contact. The magnetic 
field is studied in an arrangement of two parallel circular disks situated at 
the center of a Helmholtz configuration (Fig. |2|). This assembly of two 
identical coils with mutual distance equal to their diameter is known for the 
uniformity of the magnetic field that can be generated in it. 


In contrast to the first application where we determined A, this time the 
interest is in curl A (Eq. 4) which is obtained from linear combinations of it 
derivatives. 

There is an analogy with the derivation of stresses from displacements and sin 
the problem is essentially two-dimensional the choice for solid ring elements 
obvious. However according to SPR 933 complex stresses cannot be derived for 
these types of elements. That is why the use of three-dimensional elements may 
be considered. 

A is no longer unidirectional and Eq. (9) in fact consists of two independent 
equations, one for each of the components of X. Since no relation between thes 
components exists the equations can be solved in two subcases using the dis- 
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lacement approach. This requires that E = G and v = 0 such that the stress- 
train matrix transforms into the identity matrix multiplied by a constant, 
his stress-strain relation implies the material is anisotropic. It can be 
iPecified on MATl bulk data cards while for solid elements the material axes 
b which it refers are the axes of the basic coordinate system. 


I 


For the first subcase degrees of freedom 2 thru 6 of all grid points are con- 


3Ai 


trained. This means the a data of subcase 1 will denote E t 


zx 


= E 


3Ai 


xy ~ ^ second subcase constraining all degrees of freedom but 

3A, ^ 3A3 

T = E — ^ 


yz 


3y 


3z 
01 

and 


and 


he third at all grid points stress data a contain E 

I ^ 9A3 ^ ^ 

|Zx “ ^ ax * neglect of constant E the components of B (Eq. 4) can now 

>e obtained. The first component of rot A is found to be x given in the second 

lubcase. The second component is represented by the x data obtained from a 

ZX 

iUBCOM in which SUBCASE 1 data are subtracted from SUBCASE 2. The third compo- 

lent is the negative value of x output for the first subcase. 

! 

[nfortunately the SUBCOM results were erroneous which is presumably due to any 
if the SPR’s mentioned before. As a consequence the subtraction of results has 
;o be performed separately. Attempts to output complex stresses via OFP lead to 
'AKUNPK errors when more than one subcase is involved. Consequently x and 
are obtained from a printout of table OESl. 


I The ALTER packet required for this class of applications is a simplification 
>f the one given in Appendix I. The modifications made to it appear from the 
[ollowing: 


1 ) The FILE instruction is removed. 

1) Since the displacement approach is used the SMA2 instruction is the one 
standard in rigid format 1 . 

3) The part related to multiplication and addition of load vectors is deleted. 

+) The instructions dealing with the computation of quantities required for 
heat production calculation are no longer appropriate. 

3) The standard rigid format 1 call for SDR2 may be used but OESl has to be 
printed out. 


The resulting packet is given in Appendix II. 

I 

! 

^ CONCLUDING REMARKS 


! 

I Electromagnetic field calculations mean an area of relatively new finite 
element applications. Since the development of a production program based on 
finite element techniques has shown out to be an extensive task, those inter- 
ested in these calculations are advised to verify to which extent they can reap 
the fruits of what has already been accomplished in structural analysis. The 
convenient way of programming offered in the DMAP language and the large number 
of available features such as checkpointing, direct matrix and table input, 
parameter operations and plotting facilities makes consideration of NASTRAN 
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application worthwhile. Since many of the problems dealt with in electromagnet; 
field analysis are axisymmetric it is felt desirable to increase the 1.0 prior: 
ity status of SPR 933. 

Moreover fixes of SPR’s 458 and 483 will allow for proper execution of modules I 
ADD, ADDS and MPYAD when processing complex input matrices. Proper accounting J 
the boundary conditions and the capability of determination of quantities deri^ 
from the magnetic vector potential will be enabled by these corrections. 

It is hoped the present work demonstrates the usefulness of NASTRAN in this ar< 
of engineering and will lead to the increasing number of non-structural applicj 
tions it deserves. ; 


I 


i 


I 
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.ALTER 3 

[file 

ALTER 2?. 
[PARAM 
iALTER 35 
|SRA2 

ALTER 1*1* 
AOQ 

CHKPNT 
EO'JIV 
ALTER 47 
:SMA3 
IALTER 89 
IPARAM 
OFCOMP 

IsAVE 

CONO 
ALTER 97 
T 

I* 

PARAM 
PARAM 
PARAM 
PARAM 
CONO 
CONO 
■ AOO 
CHKPNT 
i JUMP 
LABEL 
AOO 

CHKPNT 

JUMP 


APPENDIX I 


* 4^ 

* ALTER JACKET FOR THE COMPUTATION OF THE MAGNETIC ♦ 

* VECTOR POTENTIAL IN THREE-PHASE CONOUCTING SYSTEMS * 

* * 


PG1=SAVE/PG2=SAVE S 

22 

//C»N,NOP/V,N,SKPMGG=l $ GENERATE HEAT CAPACITY MATRIX 
35 

C STM, MPT , ECPT ,GPCT,DIT/ ,MGG/V , Y, WTMASS= 1. 0/V,N, NOMGG/ 

V ,N,N03GG/C,N,-1 f 
44 

MGG,KGGX/KGG1/C,Y,0MEGA = ( 0.0,-314.1592654) I 
KGGl X NOW MODIFY TO ENABLE PROCESSING OF KGGl 
KGGl, KGG/NOGENL X EXCLUDES USE OF GENERAL ELEMENTS 
47 


89 


97 


GEI,KGG1/KGG/V,N, LUSET/7, N, NOGENL/V, N, NOSIMP S 

//C,N,NOP/V,N, counts! $ SET LOOP CONTROL COUNTER 
K LL/LLL,ULL/C,N, 1/C, M ,0/V,N,MINOIAG/V,N,O£T/V,N, POWER/ 
V,N,S1NG S DECOMPOSITION BY RBMG2 REPLACED 

MlNDIAG,OET, POWER, SING t 

HERROR3,SING $ TERMINATE COMPUTATION IF KLL IS SINGULAR 


BEGIN OF ALTER TO MPYADD LOAD 
FROM SUBCASES 1, 2 AND 3 
//C,M, ADO/ V, N,NSKIP/V ,N,N SKIP 
/ /C, N, SUB /V, N, C CUNT /V,N, COUNT 
/ /C, N, SUB /V,N,CASE3/V,N, COUNT 
//C,N ,SUB/V ,N,CASE2/V ,N ,CASE3 
LOAD3,CAS£3 $ 

L0A02,CASE2 S 
PG,/PG1 S 
PGl J 
EN03 $ 

L0AD2 $ 

PG,/PG2 S 
PG2 $ 

EN03 $ 


VECTORS obtained 

$ 

t 

$ 

S 
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APPENDIX I - Concluded 


LABEL 

AD05 

CHKPNT 

SETVAL 

SAVE 

LABEL 

RLPT 

5 

S 

EOUIV 

ALTER 100, 
SSG2 

ALTER 103, 
$ 

FBS 

CHKPNT 

SORl 

CHKPNT 

VEC 
VEC 
CHKPNT 
ADD 
HPYAO 
CHKPNT 
HPYAD 
■TPYAD 
CHKPNT 
HATFRM 
ALTER 119, 
% 

CASE 

OPO 


SOR2 

ALTER 137 
LABEL 
PPTFARM 
ENDALTER 


L0A03 f , 

PG,PG2,PG1,,/PSUM/C, Y,PHASE3= (-.5 86602 54) / 

C,Y,PHASE2=(-.5, .8660254) $ 

PSUH $ I 

//V,N , REPEAT/C, N, -1 $ 

REPEAT $ 

EN03 $ ^ 

LBL7,2 f ! 

END OF ALTER TO MPYAOD LOAD VECTORS OBTAINED FROM 
SUBCASES 1, 2 AND 3 j 

PSUH,PL/N0SET S I 

100 ^ j 

UStT, GM,YS,KFS,GO,OM,PSUM/QR,PO,PS,PL '$ 

SSG3 REPLACED BY FBS. LAST PARAMETER FOR CSP OUTPUT 

LLL,ULL,PL/ULV/C,N,1/C,N,1/C,N,2/C,N,3 % 

ULV S 

USET,PG,ULV,,YS,GO,GM,PS,KFS, KSS, QR/U GV , PGG , QG/C , N , 1 / 
C,N, STATICS $ 

PGG,UGV S 

BUILD VECTOR WITH ALL COMPONENTS EQUAL TO 1 
USET/FSET/C, N ,G/C,N, COMP/C, N,L $ 

USET/SSET/C,N ,G/C,N, L/C, N, COMP 5 
SSET, FSET S 
SSET,FSET/UNIT $ 

MGG,UGV,/TST/C,N, 0 S 
UNIT,TST t 

UNIT, TST, /ADJ1/C,N,1 S QUANTITIES ADJl AND ADJ2 ENABLE 
UGV ,TST,/ADJ2/C,N,1 S COMPUTATION OF HEAT PRODUCTION 
A0J1,ADJ2 $ 

A DJI, A0J2 , , , // $ 

119 

ADJUST SDR2 INPUT TO ALLOW FOR COMPLEX UGV OUTPUT 

CASeCC,/CAS£XX/C,N,CEIG£N/V,N, REPEAT $ 

0 YNAMICS,GPL, SIL,USET/GPLC,SILD,USETO, ,,,,,,, EQDYN/V,N, 
L USET/V,N,LUSETD/V,N,NOTFL/V, N,NODLT/V,N,NOPSOL/V ,N, 
NOFRL/V, N ,NONFLT/V,N,NOTRL/V ,N,NOEEO/C,N , 123/V,N , NOUE t 
CASEXX,CSTM,MPT,DIT,£QDYN,SILO,GPTT,EOT,BGPOT,CLAMA,QGi 
UGV ,EST,/0PG1,00G1,0UGV1,0ES1,0EF1,PUGV1/C,N,CEIG f 

HERR0R3 $ 

//C,N,-3/C,N,HSTATICS J 
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APPENDIX II 


[ 


1 $ 

$ 

'$ 

$ 

S 
S 

Ul tfp 


^ " * 
■» A|,TFP PAPKFT FOP IMF COMPI )T A T T ON OF MAGNFTTC FTFIOF » 

* 

■iMMMJ’ <M> * # -if *<;• if <S- •tMHHHnn} <}. ^HH;. <;. ,;f ^;. ^j. ^ ^ ^ ^ ^ ^ ^ ^j, ^ ^ ^ ^ ^ ^ ^ ^ ^ 


DAPAM 
A I TF^ h 
j/inn 
CHKPMT 
f^miTv 

Ialtpp ui^uy 
Ifma? 

|a| TPD PQ.BO 

nFpOMp 


FAVF 

'roMO 

;AITFP ini.lio 

I ^ 

1 

I PHKPA T 

1 sno] 


! AI TPP 1 ?1 , 1 ?1 
CMKPMJ 
! TAOPT 
A| TFP 1 ?7 


// C«N»MOP/V .N. FKPwrip = l <t GFNFPaTF MASF MATRIX 

MPG.KGPX /KGGl/C.Y.OMPGAsfO.Ot-TlO.lSQPFSA) $ 

K'PGl t- MOto PFNAMF TO F^'APl.F pdOCFFFTMG OF KGG1 

KGGl .KGG/MOGFK'L T iiGAGP OF GFNFPAI. FI f^MFMTS IF FxrLUHFn 

GFT ,KGG1 /KGG/V.N,(.i igFT/\/,N,MOGFNI /\/«N.NOfIMP J 

•^LL/I.LI «l 'L I /r »M « 1 xr *M « n/\/ » M, Mj not AG/V/ «M»OFT/V «N»POWFP/ 
\/.N.SIMG F PcplaCF P«MG? PifCOMdOF t T T ^N . KJ L POMPLFX 
MTMOI AP.nFT.POWFP.FTMGF 

HFPROPT.STMG T TC-pminATP COmpujATTOM TF KLI IF FlMGIILAp 

FFG? PFPl.ACPn RY PRS, I.AFT PAPAMFTFR FOP CFP OUTPUT 
I LL.'.'l F.Pl./IJI V/C.M. 1/C.N« l/C.M,?/r.N.T $ 

HI..V F 

IIFFT ♦PF,liLV, .YG.GO,GM.PF.KFF.KGS.OP/HGV/.PGr,,QG/\/.N, 
NFKIP/r.M.FTATTGF F 

F coMPLFx gtpffffg apf tm offi put format 

OFGl,,,,// F IF TMPRAPFP FOR fOPPFGT OFP PPOCFFFIMG 


I ARFI HFPPOPR F 

I DPTPAPM /Xr*N»-1/C*N»HFTATTCF F 

I FMHAI TFP 


I 
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